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ABSTRACT 


During  the  contract  year  1  October  1978  through  30  September 
1979  considerable  progress  was  made  toward  implementing  a  new 
Lincoln  Orbit  Determination  (LODE)  software  capability  at  MIT 
Lincoln  Laboratory  for  precise  post-mission  determination  of 
reentry  vehicle  trajectory  and  motions  during  both  exo-  and 
endo- atmospheric  phases  of  the  flight.  Early  in  the  year,  RRI 
completed  preparation  of  all  equations  for  the  Kalman-based 
forward-backward  filter-smoother  algorithm  which  provides  the 
optimal  LODE  estimation  capability.  Instructions  for  programming 
these  equations  were  prepared  in  detailed  Hierarchical-Input- 
Process-Output  (HIPO)  format  and  were  delivered  to  MIT/LL  for 
coding,  verification,  and  validation.  At  year  end,  the  equations 
have  been  coded, and  code  debugging  and  numerical  equation  veri¬ 
fication  are  in  progress,  with  completion  planned  for  calendar 
1980. 
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I.  INTRODUCTION 


This  report  sumnarizes  activities  and  accomplishments  by 
Riverside  Research  Institute  (RRI)  in  performing  metric  accuracy 
analyses  for  MIT  Lincoln  Laboratory  (MIT/LL)  from  1  October 
1978  through  30  September  1979.  The  reporting  period  covers  an 
intermediate  phase  in  the  development  of  the  new  LODE*  post¬ 
mission  trajectory  estimation  software  system  for  MIT/LL,  and 
coincides  largely  with  the  coding  implementation  at  MIT/LL  of 
equations  previously  formulated  by  RRI1.  Detailed  debugging 
and  testing  is  in  progress  at  the  close  of  the  period  while 

system  validation  and  detailed  documentation  are  scheduled  later. 

2 

for  calendar  1980. 

The  purpose  of  this  report  is,  therefore,  to  present  an 
overview  of  status  and  accomplishment  to  date.  The  archival 
function  of  documenting  detailed  technical  results,  however, 
will  be  accomplished  as  part  of  the  equation  and  software  docu¬ 
mentation  scheduled  for  calendar  year  1980. 

The  technical  goal  for  LODE  during  1979  has  remained  the 

3 

same  as  that  originally  recommended  by  RRI,  namely  implementation 
of  a  Kalman-based  forward-backward-smoother  algorithm  which  is 
optimized  for  accurate  estimation  in  the  intended  non-real-time 
post-mission  application  at  MIT/LL.  The  functional  application 
of  this  algorithm  is  shown  in  skeleton  form  in  Figure  1-1. 


♦Lincoln  Orbit  Determination 
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The  key  elements  in  the  concept  are  the  forward  and  back¬ 
ward  filters,  which  accept  measurement  data  and  produce  estimates 
of  the  trajectory  state  vector,  and  the  smoother,  vdiich  optimally 
combines  the  forward  and  backward  trajectory  estimates  in  order 
to  produce  the  single  best  estimate  of  the  trajectory  which,  at 
each  point,  is  based  on  all  of  the  measurement  data.* 

Ground  and  on-board  measurement  data  are  collected  and  pre- 
processed  to  provide  clean,  calibrated  data  for  the  filter  inputs, 
and  the  forward  and  backward  filter  runs  are  iterated  if  necessary 
to  permit  adjustment  of  the  initial  conditions,  the  parameters 
of  the  measurement  noise  models,  or  the  aerodynamic  model  which 
relates  vehicle  attitude  and  motion  to  aerodynamic  forces  and 
moments.  The  pre-processing  procedures  are  the  responsibility  of 
MIT/LL  staff  while  RR1  has  had  responsibility  for  formulating 
the  overall  algorithm  and  deriving  and  specifying  all  equations 
for  the  two  filters,  the  smoother,  the  measurement  models  for 
each  sensor  type  (i.e.,  radar,  optics.  Multistatic  Measurement  System, 
system,  on-board  gyros,  and  on-board  accelerometers),  and  the 
aerodynamic  model  (including  angle-of-attack  effects  and  vehicle 
mass  and  aerodynamic  unbalance) . 


♦Throughout,  "optimal"  and  "best"  are  used  to  mean  the  smallest- 
variance  estimate  among  all  unbiased  and  locally  linear  estimates 
of  the  trajectory  state  vector. 
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II.  STATUS  AND  ACHIEVEMENTS 


The  Metric  Accuracy  Analysis  activity  is  a  joint  effort  involving 
staff  from  RRI,  MIT/LL,  and  Ford-Aerospace,  with  areas  of 
principal  responsibility  as  shown  in  Figure  II-l.  The  status 
and  achievements  for  each  of  the  several  RRI  task  areas  will  be 
summarized  in  separate  paragraphs  which  follow.  However  the 
overall  program  status  is  as  indicated  in  Figure  II-l. 

The  preparation  and  verification  of  equations  are  essentially 
complete  while  the  completion  of  HIPO*  instructions  for  programming, 
to  convey  these  equations  to  Ford-Aerospace  programming  staff,  will 
be  completed  in  1980.  Coding,  debugging,  and  testing  are  in 
progress,  and  regular  technical  and  software  liaison  meetings 
have  taken  place  throughout  the  year  and  will  continue.  Perfor¬ 
mance  evaluation  is  the  overall  responsibility  of  MIT/LL,  with 
heavy  involvement  by  both  RRI  and  Ford-Aerospace,  and  will  be 
completed  in  calendar  1980  when  the  overall  LODE  system  (cf. 

Fig.  1-1)  has  been  debugged  with  simulated  data  and  is  running 
as  a  system.  (Subsystem  tests  of  the  filter  algorithm  thus 
far  have  been  successful  and  highly  encouraging) .  Documentation 
will  be  assembled  and  published  by  MIT/LL  in  1980  with  major 
inputs  from  RRI  as  well  as  from  Ford-Aerospace  personnel. 

The  more  detailed  description  of  status  in  RRl's  respective 
areas  is  provided  in  the  following  paragraphs,  which  in  sequence 
and  titling  are  keyed  to  the  RRI  proposal  of  activities  for 

4 

the  year  . 


♦Hierarchical  Input-Process-Output 
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FIG  H  -1  FUNCTIONAL  OVERVIEW  FOR 


METRIC  ACCURACY  ACTIVITIES 
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A.  PREPARATION  OF  EQUATIONS 

The  derivation  by  RRI  of  detailed  equations  for  the  forward 
and  backward  Kalman  filters  and  for  the  smoothing  algorithm  was 
essentially  completed  during  the  previous  contract  year.  These 
equations  provide  for  best  estimation  of  metric  trajectory  infor¬ 
mation,  ground  and  on-board  sensor  biases,  and  vehicle  aero¬ 
dynamic  parameters,  and  include  a  full  rigid-body  motion  solution 
for  a  spinning  quasi-symmetric  vehicle. 

During  the  current  contract  year  these  equations  were  extended 
to  provide  for  a  generally  asymmetric  vehicle  (that  is,  a 
vehicle  whose  aerodynamic  and  inertial  properties  are  functions 
of  roll  angle) ,  and  to  include  mass-offset  and  trim-angle  mis¬ 
alignments  between  the  principal  aerodynamic  axes  and  the 
principal  inertial  axes.  This  model  has  sufficient  generality 
to  separately  account  for  and  estimate  the  principal  physical 
origins  of  “roll-resonance"  effects.  With  this  extension  to 
a  general  aerodynamic  model,  work  on  defining  the  equations 
for  the  LODE  processor  was  complete. 

In  addition,  a  considerable  effort  during  the  year  was 
devoted  to  careful  derivations  of  alternative  forms  for  several 
critical  equations,  in  order  to  preserve  their  mathematical 
equivalence  to  the  optimal  solution  but  to  enhance  their  tractabi- 
lity  for  computer  implementation.  The  specific  critical  areas 
which  may  be  mentioned  in  this  respect  are: 

1.  propagation  of  forward  and  backward  covariance  matrices, 
denoted  P  and  Q  respectively,  for  the  filter  state 
estimates, 

2.  determination  of  their  inverses,  P-1  and  Q-*, 

3.  determination  of  the  covariance  matrix,  M,  for  the 
smoothed  estimate,  where  M  is  defined  in  terms  of 
its  inverse  by 
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4.  maintenance  of  accuracy  in  calculations  involving 
widely  disparate  magnitudes. 

From  the  outset  RRI  has  recognized  the  undesirability  of 
directly  computing  the  three  inverses*  of  very  large  matrices 
(order  100  or  greater)  which  would  be  required  for  straight¬ 
forward  implementation  of  the  smoothing  algorithm  in  its  most 
direct  form6.  instead  RRI  has  formulated  the  equations  to 
provide  for  propagation  of  the  respective  inverse  matrices 
directly,  rather  them  propagation  of  the  basic  matrix  followed 
by  direct  inversion.  The  evaluation  of  these  alternative  approaches 
has  been  accomplished  through  a  combination  of  theoretical  analysis 
(cf.  section  II-E,  "Analysis  of  Test  Cases")  and  application 
to  simulated  data  as  part  of  equation  verification  (  cf .  Section 
II-B, following,  "Verification  of  Equations").  At  the  close 
of  the  contract  period  these  evaluations  were  in  progress.  During 
preparation  of  this  report  these  evaluations  were  largely  completed. 
Only  a  definition  of  how  best  to  proceed  with  the  M  matrix  re¬ 
mained  undecided  and  was  approaching  conclusion. 

Three  RRI  contributions  deserve  separate  note.  First,  the 

exactly  optimal  equations  for  combining  forward  and  backward 

estimates  in  a  smoother  were  derived  during  the  year6  and 
7 

published  .  These  equations  extended  and  clarified  previous 
versions  existent  in  the  literature.  Second,  a  reformulation8 
of  the  basic  smoother  equation  led  to  a  form  which  reduces  by 
many  orders  of  magnitude  the  sensitivity  of  the  smoothed  estimate 
to  a  required  differencing  of  large  numbers.  Third,  an  exact  al¬ 
gorithm9  for  inverting  the  sum  of  two  matrices  has  been  derived** 


V*1,  Q_1,  and  M  *  (M-1)”1. 

**under  separate  auspices  at  RRI. 
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and  is  available  should  covariance-matrix  inversion  re-enter 
the  picture.  The  two  latter  contributions  have  far-reaching 
impact  for  large-scale  smoother  implementations  but,  to  date, 
neither  has  been  found  in  the  published  literature. 

B.  VERIFICATION  OF  EQUATIONS 

During  the  previous  contract  year  the  equations  being 
derived  by  RRI  were  verified  by  being  coded  in  a  reduced  form 
at  RRI  (the  "proto  program")  and  tested  against  simulated  data 
and  theoretical  results.  The  purpose  was  to  find  and  eliminate 
any  analytical  or  algebraic  errors  by  demonstrating  working 
subsystem  solutions  before  coding  the  full  equations  at  MIT/LL 
and  attempting  to  debug  possibly  invalid  equations.  During  the 
previous  year  the  equations  for  the  forward  filter  were  thus 
verified  and  instructions  for  programming  were  delivered  to 
MIT/LL. 

At  the  start  of  the  current  year  (November  1978)  instructions 
for  programming  the  backward  filter  and  the  smoother  were 
prepared  and  delivered  to  MIT/LL  to  be  coded  for  verification 
of  equations  in  concert  with  code  debugging.  Through  this 
method  critical  equation  areas  were  identified  and  RRI  undertook 
to  provide  either  improved  equation  approximations  or  more 
numerically  tractable  versions  of  exact  equations, as  appropriate. 

This  working  procedure  continued  iteratively  during  the  year  and 
came  to  include  code  checkout  (by  extension  of  the  proto-program) 
at  RRI  and  transfer  of  code  to  MIT/LL  via  punched  cards  and  listings. 
Indicative  of  the  analyses  and  activities  which  RRI  performed  to 
resolve  difficulties  in  equation  verification,  are  the  following: 

1.  examination  of  first,  second,  third  order  matric 
Taylor  series  for  covariance  matrix  propagation, 

2.  rearrangement  and  inclusion  of  terms  in  covariance 
equations  in  consistent  increasing  order  of  approximation, 
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3.  examination  of  alternative  and  explicitly  symmetric 
form  for  covariance  update  equation, 

4.  improved  propagation  accuracy  through  subdivision 
of  the  basic  update  time  interval  (0.1  s), 

5.  consideration  of  different  integration  schemes  for 
covariance  propagation  by  direct  integration  of  the 
covariance  differential  equation, 

6.  adoption  of  the  exact  system  transition  matrix  (compu¬ 
ted  from  the  matric  exponential  one-sided  Green's 
function)  for  propagating  covariance  matrices, 

7.  incorporation  of  mid-point  evaluation  for  transition 
matrices  to  locally  symmetrize  forward  and  backward 
propagation  equations, 

8.  incorporation  of  new  form  of  smoother  equation 
(mentioned  in  Sec.  II-A)  to  reduce  sensitivity  to  errors 
in  differences  of  large  numbers. 

This  facet  of  the  LODE  implementation  has  been  unexpectedly  drawn 
out  and  has  acted  to  delay  delivery  of  instructions  for  programming 
the  general  aerodynamic  model. 

As  of  this  writing  the  filter  equations  have  been  thoroughly 
investigateu  and  equations  for  forward  and  backward  filters  may 
be  set  down  with  same  confidence.  The  exponential  form  for  the 
transition  matrix  for  propagating  the  smoother  covariance  matrix 
has  been  verified  and  found  to  be  acceptable  when  system  noise  is 
zero.  However,  this  approximation  has  been  found  to  be  inadequate 
when  system  noise  is  present.  The  difficulty  on  the  latter  case 
was  determined  to  be  caused  by  the  fact  that  the  assumptions  used 
in  deriving  the  exponential  form  of  the  smoother  transition  matrix 
were  violated  when  system  noise  is  not  zero.  An  alternative  form 
for  propagating  M  which  does  not  rely  on  transition  matrices  is 
currently  being  investigated. 
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C.  INSTRUCTIONS  FOR  PROGRAMMING 

The  improved  programming  technique  called  HIPO*  had  previously 
been  recommended  by  RRI  and  accepted  by  MIT/LL  as  a  suitable  means 
for  transferring  a  unique  definition  of  LODE  equations  and 
algorithmic  structure  derived  at  RRI  to  MIT/LL  for  coding  and 
debugging.  Table  II-l  shows  the  schedule  of  delivery  and 
current  status.  Through  an  intensive  effort  from  5  April  1978 
in  the  previous  contract  year  to  22  March  1979  in  the  current 
year  a  set  of  HIPO  was  generated  and  delivered  by  RRI.  This 
set  (of  about  250  pages)  provides  complete  equations  for  LODE 
with  only  three  exceptions. 

First. equations  for  the  general  aerodynamic  model  have 
been. derived  at  RRI  but  reduction  to  HIPO  form  remains  to  be 
done.  Second,  nose-tip  range  and  range  rate  extensions  to  the 
measurement  model  remain  to  be  derived.  Similarly,  and  third, 
the  estimation  equations  for  (small)  timing  bias  are  planned 
for  definition  and  phase-in  in  1980. 

In  addition  RRI  has  provided  improved  and  revised  covariance 
propagation  and  smoother  equations  to  MIT/LL  in  the  form  of 
card  decks  and  listings  to  reflect  improvements  discussed  above 
in  Section  II-B.  It  is  expected  that  replacement  HIPO  sheets 
will  be  provided  to  document  these  and  any  further  equation  or 
algorithm  changes  which  may  be  necessary.  Also,  temporary 
resource  constraints  at  RRI  are  requiring  the  delay  of  further 
effort  on  HIPO  because  of  the  higher  priority  which  has  been  assigned 
to  investigation  of  the  smoother  algorithm. 

D.  ANALYSIS  OF  TEST  CASES 

Several  major  different  forms  of  test  cases  have  been  used 
to  debug,  verify,  and  checkout  coded  equations. 

♦Hierarchical  Input-Process-Output 
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At  RRI,  simulated  data  corresponding  to  point-mass,  zero- 
drag  motion  in  a  central  inverse-square  force  field  was  used 
extensively  early  in  the  year  to  analyze  filter  performance  and, 
as  a  result,  to  find  and  eliminate  anomalies.  The  explanation 
of  initially  questionable  trends  in  (forward)  estimator 
covariance  was  eventually  found  in  aspect  angle  effects  at  low 
altitudes,  and  in  the  initial  effects  of  (zero)  drag  uncertainty 
at  higher  altitudes.  Also,  estimation  residuals  were  materially 
reduced  (and  brought  into  consistency  with  optimal  covariances) 
when  effects  due  to  non-linearity  in  the  measurement  model  were 
identified  and  the  sequence  of  processing  radar  measurements 
was  changed  frcm  R,  E,  A  to  A,  E,  R. 

Also,  a  case  corresponding  to  one-dimensional  particle 
motion  along  a  straight  line,  with  constant  drag  parameter,  was 
formulated  and  solved  to  achieve  analytic  expressions  for  the 
covariance  matrix  elements.  These  results  were  used,  in  con¬ 
junction  with  the  previous  simulation,  to  achieve  considerable 
confidence  in  the  filter  covariances  and  in  very  large  allied 
sections  of  the  equations  and  the  system  and  measurements  models. 

Third,  simulated  input  data  corresponding  to  a  known 
ballistic  coefficient  profile  are  being  used  in  analyzing 
smoother  algorithms,  where  the  assumption  of  non-zero  process 
noise  is  vital. 

Finally,  detailed  hand  calculations  have  been  performed 
by  RRI  down  to  the  individual  subroutine  level  for  the  HIPO 
which  has  been  delivered  to  MIT/LL.  These  calculations  provide 
a  valuable  resource  which  can  be  useful  in  establishing  confidence 
in  detailed  LODE  calculations  after  subroutine-by-subroutine 
comparison  between  MIT/LL  and  RRI  results  . 

E.  TECHNICAL  AND  SOFTWARE  LIAISON 

During  the  contract  year  liaison  among  RRI,  MIT/LL,  and 
Ford- Aerospace  was  achieved  through  periodic  (approximately 
monthly)  project  review  meetings,  supplemented  by  smaller  working 
meetings  at  more  frequent  intervals  as  necessary.  Typically, 
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the  delivery  of  HIPO  was  accompanied  with  a  session  devoted 
to  a  detailed  walk-through  by  RRI  of  the  new  HIPO,  plus  the 
resolution  of  any  questions  for  the  previous  HIPO.  Topics  which 
have  formed  the  subjects  of  specific  meetings  include: 

1.  discussion  of  HIPO  for  processing  Multistatic 
Measurement  System  data, 

2.  RRI  recommendations  for  creation  of  source  of  six 
degree-of-freedcm  (6  DOF)  simulated  data  at  MIT/LL, 

3.  use  of  already  delivered  rigid  body  HIPO  as  a 
kernel  for  the  MIT/LL  six  degree-of-freedom  simulator, 

4.  explanation  of  equations  in  card  decks  providing 
new  covariance  propagation  equations  to  MIT/LL. 

F.  LODE  PERFORMANCE  EVALUATION 

Evaluation  of  the  performance  of  the  LODE  system  with  respect 
to  such  factors  as  accuracy,  convergence,  input  parameter 
selection,  aerodynamic  modeling,  measurement  noise  and  sensor 
selection,  assumed  process  noise  parameters,  and  run  iteration  pro¬ 
cedure  is  contingent  upon  running  the  system  as  a  whole.  However,  that 
was  not  achieved  for  an  entire  trajectory  during  the  contract 
year. 

Nevertheless,  examination  of  specific  subsystem  behavior  (e.g., 
forward  filter)  was  possible,  as  well  as  observation  of  entire 
system  behavior  over  restricted  regimes  (e.g.,  exo-atmospheric 
with  zero  process  noise) ,  and  qualitatively  encouraging  results 
and  expected  behavior  observed.  The  forward  filter,  for  example, 
provides  estimates  which  are  consistent  with  the  input  data  and 
shows  covariances  which  decrease  in  close  agreement  with  theory. 
Smoother  estimates  (with  zero  process  noise)  are  compatible  with 
the  input  estimates  from  the  two  filters,  and  have  covariances 
which  throughout  the  trajectory  are  smaller  than  both  filter 
covariances  (as  expected) . 
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Further  quantitative  performance  evaluation  has,  however, 
been  given  lower  priority  until  a  properly  working  smoother 
can  be  certified. 
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III.  STATUS  SUMMARY 

During  the  contract  year  rri  completed  the  preparation  of 
equations  for  the  full  rigid-body  general-aerodynamics  LODE 
capability,  and  delivered  complete  instructions  for  programming 
the  full  capability  with  quasi-symmetric  vehicle  aerodynamics. 
Augmentation  of  these  instructions  for  programming,  to  include 
the  general  aerodynamics  package, is  planned  for  calendar  year  1980. 
Both  point-mass  and  rigid-body  codes  have  been  produced  at  MIT/LL 
from  the  RRI  instructions  for  programming  and  intensive  code  and 
equation  debugging  are  now  in  progress  for  the  simpler  point- 
mass  code.  The  tasks  which  now  remain  to  be  accomplished  to 
achieve  the  6  deg ree-of- freedom  LODE  capability,  in  addition  to 
equation  and  code  documentation,  are  shown  in  Figure  III-1. 

During  preparation  of  this  report,  a  three  degree-of-freedom 
(3  DOF)  LODE  capability  has  been  debugged  to  the  point  of  generating 
complete  forward  and  backward  filter  estimates  and  computing 
optimal  smoothed  estimates  for  metric  data  and  ballistic  coefficient 
throughout  reentry.  This  marks  completion  of  a  major  milestone, 
and  further  testing  with  lift  forces  and  bias  estimation  may 
now  proceed  with  the  code  already  prepared.  AI90,  as  indicated 
in  Figure  II-l,  simplified  aerodynamic  coefficient  driver  equa¬ 
tions10  have  been  delivered  to  MIT/LL  to  permit  simulation  of 
rigid-body  aerodynamic  motions  and  trajectories.  The  RRI  task 
for  preparing  additional  instructions  for  programming  the  general 
aerodynamics  package  has  been  in  abeyance,  as  full  RRI  resources 
have  been  devoted  to  achieving  the  working  point-mass  solution 
and  code  already  mentioned. 
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FIGlIH-1  INTEGRATION  OF  RECENT  ACTIVITIES 
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